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Abstract 

We argue that the Mellin-Barnes representations of Feynman diagrams can be 
used for obtaining linear systems of homogeneous differential equations for the orig- 
inal Feynman diagrams with arbitrary powers of propagators without recourse to 
the integration-by-parts technique. These systems of differential equations can be 
used (i) for the differential reductions to sets of basic functions and (ii) for counting 
the numbers of master integrals. 
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1 Introduction 



The importance of the differential-equation approach to the description of the an- 
alytical properties of Feynman diagrams has been recognized a long time ago pQ. 
Within dimension regularization [2] , the differential-equation technique [3] , based on 
the integration-by-part (IBP) relations [3], has been one of the most popular tools 
for the analytical evaluation of Feynman diagrams during the last decade [5] . Here 
we argue that linear systems of homogeneous differential equations may be derived 
for Feynman diagrams starting from their Mellin-Barnes representations without 
resorting to IBP relations [6]. 

Our staring point is the multiple Mellin-Barnes representation of Feynman dia- 
grams [7 8 9j, which may be written in the following form: 



/+ioo 
dt^z 1 
-ioo 



Const ( l °° TT dt r( S^i Agjti + Bg) E ; agti m 
Const/ H dtc nZUC bj t J+ D b ) Z * ' (1) 



100 a,b,c,r 



where Zk are ratios of Mandelstam variables and A,B,C,D,a are matrices and 
vectors depending linearly on the dimension n of space-time and the powers of the 
propagators. An important property of Feynman diagrams is that the matrices A 
and C only include integers. Let us define the polynomials Pi and Qi as 

Pi® = 0+ el) (2) 

Qi® Hi) ' 

where e~l is a unit vector with nonzero element at the i-th place. Then, the integral 
in Eq. ([1]) satisfies the following linear system of homogeneous differential equations: 

Q i (% ^ e -^A,B;C,D-^ = P i (r)\ t ^ 9 ^A,B-,C,D-,g) , (3) 

where Qi = Zid/dzi, we assume that, in Eq. ([1]), the number of variables Zi is equal to 
number of integration variables ij, and, for simplicity, we put aij = 5ij. Moreover, 
the function $ defined in Eq. ([T|) satisfies the differential contiguous relations: 

$(A,B+e* a ;C,D;z)= (^AaA+B^j ^(A, B; C, D; z) , 



$(A,B;C,D-e b ;z) = [ Y^C bj 6 j+ D b ] $(A, B;C, D; z) . (4) 



This system of differential equations may be analyzed using the Grobner basis tech- 
nique |10| . Specifically, (i) the holonomic rank r of the system in Eq. ([3]) can be 
evaluated^] and (ii) starting from the direct differential operators in Eq. (j4|) and 



1 This can also be done using a finite number of prolongations and converting the original system into 
a Pfafnan system [TT] , 
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the system of differential equations in Eq. (J3j) > the inverse differential operators can 
be constructed [12]. By the action of such differential operators on the function <£, 
the value of any parameter Bj^D^ may be shifted by an arbitrary integer. This 
algorithm is working for holonomic functions if the number of linearly independent 
solutions is finite. As was shown in Ref. [13], within analytical regularization |14j . 
Feynman diagrams satisfy holonomic systems of linear differential equations under 
the condition that all particles have different masses. In fact, this statement is the 
basis of the algorithm proposed in Ref. [15]. To our understanding, a rigorous proof 
for the case of on-shell diagrams or diagrams with zero internal masses does not 
yet exist, albeit very interesting work has been done in this direction [16J. For our 
analysis, it is sufficient to assume that there is a set of independent variables z% for 
each Mellin-Barnes integral! in Eq. ([1]). Moreover, it has been shown |18] recently 
that Mellin-Barnes integrals satisfy systems of differential equations corresponding 
to Gelfand-Kapranov-Zelevisky hypergeometric equations [19] . Another necessary 
condition is that two contours differing by a translation by one unit along the real 
axis are equivalent. From the analysis performed in Refs. [8; 20], we expect that this 
statement is valid for all Feynman diagrams before constructing e expansions. 

The aim of the present paper is to illustrate this approach and to outline how 
it can be used to count master integrals, considering as examples one- fold Mellin- 
Barnes integrals and the corresponding Feynman diagrams. In fact, the evaluation of 
master integrals is often the bottleneck of multi-loop calculations of both Feynman 
diagrams and scattering amplitudes, especially if several different mass scales are 
involved, and any opportunity to reduce their number below the achievements of 
the standard techniques of solving IBP relations, as implemented in various publicly 
available computer codes, should be highly welcome. The technique advocated here 
may also allow one to gain deeper insights into the mathematical structures of 
multi-scale Feynman diagrams. 



2 One-fold Mellin-Barnes integral 

Let us consider the function 



J 1+ iR \\ k ._ 



n^ r(^ + t)njur(C7 3 --t) 
r(fl fc + t)niUr(A-t) 



B; C, D; z; r) = C* / dt x "~" J 7 ' -T(-t)z t+r , (5) 



where C$ is some z-independent constant depending on ratios of T functions with 
arguments being linear combinations of powers of propagators and the space-time 
dimension n. In the remainder of this letter, we assume that r = 0. In order to 
restore a nonzero value of r, it is sufficient to substitute — >• 6 — r. Let us assume 
that the differences between any two parameters Ai, Bj,Ck, Di are not integers. 



2 The application of the summation technique to Mellin-Barnes integrals with Zi = 1 was discussed in 
Ref. 02]. 
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Then, this function satisfies the following homogeneous differential equation: 



R L 

{-i) L+1 e Y[(e+Bi-i) Hie-CjMl, b- c, d- z ) 

i=l j=l 
K J 

= (-l) J zl[(e+A i )l[(e-D j + l)^(A,B;C,D;z) , (6) 

i=l j=l 

where 9 = zd/dz. Let us consider the case of the non-confluent function, for which 
the orders of the differential equations on the l.h.s. and r.h.s. of Eq. are equal 
to each other, viz. 

K + J = l + L + R = p, (7) 

so that the function $ satisfies a differential equation of order p. In this case, there 
are p linearly independent solutions of the differential equation. In accordance with 
Takayama's algorithm |12j . the differential operators inverse to the operators defined 
by Eq. @ can be constructed, and the result of the differential reduction applied 
to the function $ has the following form: 

p-i 

P ^l+m 1: B + m 2 ;d + m 3 ,D + m A -,z) = s ^R i 6 i <$>{A,B-C,D-z) , (8) 

i=0 

where Po,Ri are some polynomials and mj are sets of integers. 
Theorem 1: 

Any Feynman diagram associated with the function <3? defined by the Mellin-Barnes 
integral in Eq. ([5]), under the conditions that all parameters as well the differences 
between any two parameters are not integer and Eq. ([7]) is valid, has p master in- 
tegrals (including all integrals following from the original one by contracting one or 
more lines), where p is defined by Eq. (|7|). 

This theorem follows from Takayama's analysis of the differential equations in 
Eq. ([6]) with irreducible monodromy groups. 

Under the conditions of Theorem 1 and the additional condition C{ ^ Cj, it is 
possible to close the contour of integration [8] and to write the function <I> in terms 
of linear combinations of hyper geometric functions whose series representations are 
well defined in the vicinity of the point z = 0: 



z ~\(A*nn^ n£ir(^)nt=i r (^) p (A,l-D 

— — <J>( A, B: C, D: z) = — B — 7 \ -> -> 



z c m n^=i r (Aj+c m ) Y[j = i-j^ m F (Cj -C„ 
' ^ " nf=ir(5 fe +c m )niLir(A-c m ) 



m=l 

A + C m , 1 + C m — D 
B-\-C m , l-\-C m —C 



x pFp-i 



{-l)^ J z ] , (9) 



where C denotes the set of parameters C excluding C„ 
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Corollary 1: 

The number of nontrivial master integrals of a Feynman diagram satisfying the 
conditions of Theorem 1 and the additional condition that C, 7^ Cj is equal to the 
number of basic functions for any hyper geometric function on the r.h.s. of Eq. ([9]). 
Comment 1: 

The application of the relation 

T(nx) = n nx ~ 1 ' 2 {2T:) l ^T{x)T (x + ^j---v(x+ , (10) 

where n is integer and x = t + p/q, allows us to express any one- fold Mellin-Barnes 
integral as an integral of the form defined by Eq. ([1]). 

A special consideration is necessary when some parameters or some differences 
between parameters are integer. In this case, the final expression for the differential 
reduction of Horn-type hyper geometric functions has a simpler form (for details, see 
Ref. [21). To evaluate the dimension of the solution space of the solutions of the 
differential equation in Eq. ([6]), a classical technique can be applied [22]. However, 
for practical applications to Feynman diagrams, we only need two particular cases, 
namely case (i) when the differences between upper (Ai,Dj) and lower Cz) 
parameters are positive integers and case (ii) when some of the parameters are 
positive integers. 

In case (i), the original Mellin-Barnes integral in Eq. ([5]) may be simplified as 

/ ^ T ^T(A™t)^ F ® = ( A + m - 1+e )( A + m - 2+e ) ■ ■ ■ ( A+e ) J dtz'Fit) .(11) 

In this case, the order of the differential equation in Eq. © is p — 1, and there 
are p — 1 nontrivial master- integrals. After the application of step- up/step-down 
operators to the r.h.s. of Eq. ©, we obtain the hyper geometric functions P _i-F p _2, 
and Corollary 1 is valid. 

In case (ii), the differential equation in Eq. ([6]) has a common differential factor 9, 
since z (0+1) f(z) = (zf(z)) . In this case, the result of the differential reduction 
of Eq. ([8]) has the following form: 

p-2 

P ^A + mi,B + m 2 \C + m 3 ,D + rh4;z) = ^2R i 6 i $(A,B;C,D; z)+R p _ 1 (z) , 

i=0 

(12) 

where Pq, Ri are some polynomials and m; are sets of integers. The same is true for 
the r.h.s. of Eq. Q, so that Corollary 1 is again valid. 
Corollary 2: 

The number of nontrivial master-integrals of the Feynman diagram associated with 
the function defined by Eq. ([5]), under the condition that Cj 7^ Cj, is equal to 
number of basic functions for any hypergeometric function on the r.h.s. of Eq. ([9]). 
Comment 2: 

The differential equation in Eq. ([6]) can be factorized due to the relation 

z{9 + l+a)f{z) = (6+a)(zf(z)) , (13) 
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where a is an arbitrary parameter. This equation corresponds to the factor T(a — 
t)T(l — a + t) in the numerator or denominator of the Merlin-Barnes integral in 
Eq. ©. 
Comment 3: 

Also Tarasov [15] proposed that, in the framework of Ref. [13J, there is a one- 
to-one correspondence between the number of master integrals obtained from the 
IBP relations and the dimension of the solution space of a corresponding system of 
differential equations. 

In this way, we showed that our conjecture presented in Ref. [21] is correct. 
Below, we present its multivariable generalization: 
Proposition 1: 

When a multivariable Mellin-Barnes integral can be presented as a linear combina- 
tion of multivariable Horn-type hypergeometric functions with rational coefficients^] 
about some points zi = zf, the holonomic rank of the corresponding system of linear 
differential equations is equal to the holonomic rank of any hypergeometric function 
in its hypergeometric representation. 

The proof of this proposition is based on the same technology, namely the com- 
parison of the holonomic rank of the system of differential equations in Eq. ([3]) 
with the holonomic rank of each term of its hypergeometric representation. A two- 
variable example of this statement was presented in Ref. [23J . 
Conjecture 1: 

Any polynomial (rational) solution of a multivariable linear system of differential 
equations related to a Feynman diagram can be written as a product of one-loop 
bubble integrals and massless propagator or vertex integrals. 

3 Feynman diagrams 

Several examples of Feynman diagrams corresponding to the function <3? in Eq. ([5]) 
were presented in Refs. [21|23|24] , The results of the analysis performed in 
Refs. [21)23)24] are in agreement with Corollary 2 of the present paper. Neverthe- 
less, for the illustration of the advocated technique, we consider here the diagrams 
depicted in Fig. [TJ In all examples, we put rj = 0, since non-zero values of rj do 
not affect the order of the differential equations and may be easily restored by the 
redefinitions 0i — > 0i — rj. 

3.1 One- loop vertex QEDi 

Let us consider the one-loop QED vertex diagram with one massive internal line 
and two external lines being on mass shell, which is labeled QEDi in Fig. [TJ It is 

3 We called such a representation hypergeometric. 
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Figure 1: The diagrams on the r.h.s. emerge from the original one on the l.h.s. with 
arbitrary positive powers of propagators by the systematic contraction of one line. 



given by 

Cqed^o-i,^, jz) = j 



d n k 



[{k - pi) 2 ] CT i [(k + p 2 ) 2 }^ [k 2 - m 2 ]i3 
^-"W^-m 2 )^ 12 ^' 3 



pf=p 2 .=m 2 



rioo / 
™ J -too V 



U2 \ * 



(14) 



n-t) 



r (n-£7i2-j 3 )r (j 3 )T (ai)T (a 2 ) 2m 
x r( CTl +t)r(a 2 +t)r (a 12 +i 3 -|+*) r (n-j 3 -2a 12 -2t) , 

where o\ 2 = a\ + a 2 and k 2 = (pi — p 2 ) 2 . This diagram corresponds to a <E> function 
with six parameters, 

CQEDi (0-1,0-2,^3) = $ QEDi (cti, a 2 ,a 12 +j 3 - §; -; - a 12 , n ~ J 2 3+1 - ai 2 ,0; , 

(15) 

and satisfies a third-order differential equation, 

[0(0-Ci)(0-C 2 ) + zie+A^e+A^e+Az)} <D QEDl (Ai, A 2 , A 3 ; -; C U C 2 , 0; -; z) 
= . (16) 



With the help of differential operators, 

$QE Dl (^ + 1) = (0 + Ai) $ QEDl (A 



(17) 
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it is easy to reduce the value of A\ (or A2) to unity, and the one of A3 to 1 — C\ 
or 1 — C2, depending on whether j'3 is even or odd. In this case, there is a double 
factorization of the differential equation in Eq. ()16p due to Eq. (|13p , 



9(6 -d) [(0-C 2 ) + z(9+A 2 )] $ QEDi (1,A 2 , 1 - Ci; -; Ci, C 2 , 0; -; z) = . (18) 

As a consequence, there are two polynomial solutions. We conclude from our analy- 
sis that there are three master integrals corresponding to the order of the differential 
equation in Eq. (j 1 6 P : one is a nontrivial function and two are polynomials. 



3.2 One-loop vertex QED 2 

Let us now consider the one-loop QED vertex diagram with two massive internal 
lines and two external lines being on mass shell, which is labeled QED2 in Fig. [TJ 
It is given by 

d n k 



CQED 2 0'l, j2,cr) = J 



[(k - pi) 2 - m 2 ]-? 1 [(k + P2) 2 - m 2 p' 2 (k 2 ) c 



p\=p%-- 



i l-n^n/2(_ m 2)2-ii2^ r ( n _ Jl2 _ 2cJ ) 1 y*» 



r(ii) r(j 2 )r(n-ii2 -a) 2m J_ iOQ ~" V m 2 

r(-t)r(i! + t)r(j 2 + t)T (ii 2 +(T-§ +t) 



(19) 



r(ji2+2t) 

This diagram corresponds to a $ function with six parameters, 

C Q ED 2 0'i ) i2 ) ff) = *QED 2 (ii ) j2,Ji2+^-f; i 2 1 , ii | tl ,0;-;-;^ . (20) 



Repeatedly applying differential operators, <&QED 2 (-4i + 1) = (0 + A4) $qed 2 (^i) , 
and <I>qed 2 (B; - 1) = (0 + -B»-l) < &qed 2 (Bi) , we reduce Ai, A 2 , and i?i to unity, 
so that the six-parameter $ function reduces to a four-parameter one, 

$qed 2 (1, 1, A 1, 5; 0; -; z) = $ QE d 2 (1, A; B; 0; -; z) , (21) 

which satisfies a second-order differential equation, 

6 l(6+B - 1) + z(0+A)] $qed 2 (1, A B; 0; -; z) = . (22) 

Consequently, there are two master integrals, one of which is polynomial. This 
result does not depend on whether a is integer or not. 

3.3 Two-loop sunset diagram J i2 

Let us now consider the two- loop sunset diagram from Ref. [25], which is given by 

cT(kik 2 ) 



J i 2 (<7,a,/3) = 



Pi - p) 2 ] a [kf - M 2 ] a [(h - k 2 ) 2 - m 2 f 



[ i l-n 7r n/2]2(_ m 2jn-o- ^ r (n_ CT ) j poo / M 2 \ < 



r (a) r (a) T(/3) 2vri,/_ ioo V™ 2 

r (a+a-%+t) r (a+/3+o--n+i) r(-i)r (|-a-t) T (2n-a-2(x-2p-2t) 
r(n-fj-Q-t)r(^-cj-a-/3-t) 
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This integral corresponds to a f function with eight parameters and satisfies a 
fourth-order differential equation, 



(fl-f+a) (0-n+q +CJ+/3) [O-n+^+a+0) $j 012 



(24) 



: |0-~+<r+a) (0-n+a+/3+cx) (6>-n+cr+a+l) ^0- y +cr+a+/3+lj $j 012 . 

With the help of step- up/step-down operators, this equation can be written as 
follows: 



n 



■ + h 



-n+h) 



1 

-n+-+h 



3n 



-h 



. (25) 



Consequently, there are four master integrals, two of which are polynomial. All 
topologically possible integrals for this diagram are depicted in Fig. [TJ The right- 
most one is polynomial in this approach. In order to have four master integrals, it 
is necessary that the diagram with the original topology has three master integrals, 
which we may take to have the propagator powers (1,1,1), (1,1,2), and (1,2,1). 
We note that this is in accordance with the results of Ref. [26]. However, one of 
these three master integrals should be polynomial. This polynomial solution was 
first derived in Ref. [27J. An alternative derivation has recently be presented in 
Ref. EHl. 



3.4 Three-loop vertex diagram F 

Let us consider the vertex diagram denoted as F 3 in Fig. [U which is given by 

F(<Ti,<T2,<T3,ai,Q 2 ) (26) 

d n (hk 2 k 3 ) 



[{k-pi) 2 }^[{k+p2) 2 }° 2 [kl-m 2 } a ^ [k%-m 2 ] a2 [(k 1 + k 2 + k 2 )]rt 



p\=pI=o 



where we have introduced a "dressed" massless propagator, as in Eq. (45) of Ref. [2T] , 
Instead of three massless lines, we consider q%, q 2 , and q 3 massless lines correspond- 
ing to propagators with powers oi, cr 2 , and o 3 , respectively. In the present case, this 
is equivalent to the redefinition aj — > Gj — §(9/ — 1) • The Mellin-Barnes integral 
for the vertex diagram of Eq. (|26[) has the following form: 

= const x — r°° dt ( *\ , n<*+w<**n*»-\+m\+*) 



(27) 



2vri J-ioo V m 2 J T (a\2 + 2t) 

T (al-f q 3 -t) r (^ qi3 -a 13 +t) r (2<fc 3 -*23+t) r (a 123 -| (9123-1)-*) 



r (f (q 3 +i)-(T3+t) r (§9123-^123+0 

where Const is some constant irrelevant for the differential equation. 

This diagram corresponds to a $ function with twelve parameters. The first 
simplification arises from the step 

r(Qi+t)r(a 2 +t) r(/i+t) 



r(« 12 +2t) r(i+ / 2 +t) ' 



(28) 
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where I a are integers. The further procedure strongly depends on whether the values 
of qi and qi are even or odd. For simplicity, we put q\ = q2 = q 3 = 1 and denote 
this integral as F 3 . Then the following step is valid: 



r (§913-^13 + *) r (§</23-<?23 + t) 



T (%{q 3 + l)-a 3 +t) 



T(n+I 3 + t) , 



(29) 



51=92=53=1 



so that F 3 satisfies a homogeneous differential equation of fourth order that may be 
written in the following form: 



6+-n+h 



1 n T 



»+n+/ 4 )$ F3 



n 



z (e+h) [9--+i 7 



n 



) + -+h)(0+n+I 8 )<£> F3 



(30) 



Applying step- up/step-down operators, this equation can be converted to the form 



0+ 2J W+n) 



9+-n + I 2 



>+/3 



n 

h/4 

2 



<f Fs = . (31) 



Consequently, for the vertex diagram F 3 there are four master integrals, two of 
which are polynomial. All topologically possible integrals for this case are shown in 
Fig. [TJ In order to have four master integrals, it is necessary that the diagram with 
the original topology has two master integrals. 



4 Discussion and conclusions 

We proposed a novel way of finding linear systems of homogeneous differential equa- 
tions for Feynman diagrams with arbitrary powers of propagators. It is based on the 
Mellin-Barnes representation and does not rely on the IBP relations [1] . Systems of 
equations such as Eq. ([3]) are the left ideals in the Weyl algebra of linear differential 
operators with polynomial coefficients. Exploiting the Grobner basis technique |10j . 
the original diagrams may be explicitly reduced to sets of basis functions, and their 
holonomic ranks, i.e. the numbers of their linearly independent solutions, may be 
evaluated. 

For the one- variable case, we presented a very simple algorithm for counting poly- 
nomial (rational) solutions of differential equations arising from the Mellin-Barnes 
representations associated with Feynman diagrams. This reduces the problem to 
the one of factorizing the homogeneous differential equation over the field of poly- 
nomials. We presented a few simple examples to illustrate our technique. The 
generalization of this algorithm to the multivariable case is nontrivial. 

With the help of the new technology presented here, we proved the conjecture 
formulated in Ref. |21] regarding the counting of the numbers of master integrals via 
hyper geometric representations. This result may be useful for searching polynomial 
(rational) solutions of multivariable Feynman diagrams. 

We suggest that any polynomial (rational) solution corresponds to a product 
of one-loop bubbles and massless single-scale diagrams with coefficients that are 
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products of Gamma functions (see Conjecture 1). Even in the one- variable case, 
such a correspondence does not emerge from the application of standard IBP rela- 
tions, as was pointed out in Refs. |27|28j for the case of the two-loop sunset diagram 
^012 • With help of the technique presented here, all such algebraic relations between 
master integrals of the type studied in Ref. [21)23124] may be easily derived. 

Hypergeometric functions provide us with a remarkable tool for deepening our 
understanding of the mathematical structures underlying Feynman diagrams, and 
the present analysis allows us to draw the following picture. From the fact that 
any Feynman diagram with arbitrary powers of propagators is reducible to a set 
of master integrals, including bubble diagrams and massless propagators, it follows 
that a given Feynman diagram corresponds to a special function with a reducible 
monodromy group (see also Ref. |29|30j ). The dimension of the irreducible part of 
the monodromy group, which is equal to the dimension of the solution space of the 
Pfaff system of differential equations, is equal to the number of master integrals 
generated via IBP relations, provided the full set of the latter is exploited. It is 
interesting to note that the simplest way of avoiding the reducibility of the mon- 
odromy group is to introduce different non-integer parameters for each propagator. 
This may be regarded as a generalization of the analytical regularization. 

The proposed method to derive differential equations from the Mellin-Barnes 
representations of Feynman diagrams with non-unit values of z% can be directly 
applied to study massless propagator diagrams at higher loop orders. In this case, 
we have to tackle with multiple (threefold and higher) Mellin-Barnes integrals. We 
shall return to this issue in a future publication. 
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